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Abstract 

Recent experimental results for the magnetic hyperfine field B^ at the 
nuclei of s-p impurities such as 119 Sn in intermetallic Laves phases RM 2 (R = 
Gd, Tb, Dy, Ho, Er; M = Fe, Co) and 111 Cd in RCo 2 , the impurity occupying 
a R site indicate that the ratio -Bhf//-*3d exhibits different behavior when one 
goes from _RFe2 to RCo 2 . In this work, we calculate these local moments and 
the magnetic hyperfine fields. In our model, B^ has two contributions: one 
arising from the R ions, and the other arising from magnetic 3<i-elements; 
these separate contributions allow the identification of the origin of different 
behavior of the ratio mentioned above. For m Cd in RCo 2 we present also the 
contributions for B^ in the light rare earth Pr, Nd, Pm, Sm compounds. For 
the sake of comparison we apply also the model to 111 Cd diluted in i?Ni2. 
Our self-consistent magnetic hyperfine field results are in good agreement 
with those recent experimental data. 

Keywords: Magnetic hyperfine field; Laves compound; formation of local 
magnetic moment 



1. The model 

The local magnetic moment at a Sn impurity diluted in _RFe2 and RCo 2 or 
Cd in RCo 2 , has the following contributions: one from the R ions, consisting 
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of a polarization of the s-p impurity level by 4f and 5d rare-earth electrons 
and from a potential (V) due to the presence of the impurity. Also included 
in V, is the nearest-neighbor hopping change due to the break of translational 
invariance by the impurity The other one, from the transition metal, is a 
magnetic field produced by the neighbors of the impurity, which occupies a 
rare earth site (see Ref.( [1,]) and references therein). 

The Hamiltonian to describe the formation of the local s-p magnetic mo- 
ment and hyperfine field is 

U = H R + V + H M . (1) 

In Eq. ©, 

H R = J2 £ « c l c ™ + Yl **j4t c j<t» (2) 

defines an effective pure rare earth host which consists in a conduction s-p 
band polarized by the 4f and 5d electrons. In Eq. (j2J), is the center of the 
s-p energy band, now depending on the spin a orientation, c\ a (q ct ) is the 
creation (annihilation) operator of conduction electrons at site i with spin a 
and tij is the electron hopping energy between neighboring i and j sites. The 
second term of Eq. (jTJ) is the potential due to the presence of the impurity 
at site i — 0. 

V = Yl V o° c o* c Ov + ^Y tm ( c o° c ^ + c t c o<x) 7 ( 3 ) 

where V 0u = (e l 0a — e^) is a spin dependent local term, e\ a being the s-p 
impurity state energy level. The parameter r takes into account the change in 
the hopping energy associated with the presence of the impurity (H0,0],t = 
meaning no disorder in the hopping. 
The last term of Eq. (jTJ, 

H M = -J2^ Sd (S M )4^, (4) 

is the interaction energy between the magnetic field from the M ions and 
the impurity spin. J sd is an exchange parameter and (S* M ) is the average 
magnetic moment at M sites. 

Let us first solve the problem with Hm = 0. Using Dyson equation, the 
exact Green's functions 'g i j a {z) due to the charge perturbation at the origin, 
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Here 9ij a {z) is the Green's functions for the pure host with z = e + iO, and 
a = t + 1. The effective potential V° s (z) is given by: 

V: s (z)=V 0(T + (a 2 -l)(z-e h a )., (6) 

For the full Hamiltonian ([!]), the perturbed Green's functions Gij a (z) is 
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where V£ n = —aZ nn J sd (S M ) is the magnetic coupling between the impurity 
and the M sites and Z nn is the number of the nearest neighbor M ions sur- 
rounding a s-p (Sn or Cd) impurity. At this point we adopt the following 
approximation: one consider only the contribution of nearest neighbor M 
sites; and in this sense we have a cluster-like approach to calculate the con- 
tribution to the magnetic moment arising from M ions. We assume that V° n 
is small compared to Vo a thus justifying the Born approximation, T{[ « V° n . 
The local Green's functions Goo a (z) at the impurity site is then 
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Note that in the local Green's functions Goo a (z), the first term is due to the 
rare-earth ions whereas the second comes from the polarization produced by 
the M ions. The local potential Vo a is self consistently determined using 
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Figure 1: Calculated local density of states Pq(s) (a) at the Cd and (b) at the Sn impurity 
site both in GdCo2. Due to the splitting and the band deformation caused by the scattering 
produced by the charge difference between the impurity and host, these local DOS's will 
generate local magnetic moments in the indicated directions (see text). 
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the Friedel screening condition (see for instance Ref. |3J) for the total charge 
difference AZ between impurity and rare earth atoms. 

Assuming that the screening of this charge difference is made by the s-p 
band then AZ = AZ^ + AZ^, where AZ a is 

AZ a = -ilmln [a 2 - g 00a (s F ) V&(e F )] , (10) 

7T 

where e F is the Fermi energy level. The local s-p density of states per spin 
direction at the impurity site are calculated by p a (s) = (— ImG r ooo-(^)- 
The local s-p electron occupation number, n 0a , is obtained by integrating 
the corresponding local density of states up to the Fermi level £p- The total 
magnetic moment (mo) at a s-p impurity, given by m~o = not — n oi is 

m = m£ + m nd , (H) 

where 

rh R = 

° ^^i-oo"" a* - g 00a (z)V£{z) 

is the contribution from rare-earth ions and 

1 /■«* a 2 Z nn J sd (S M ) 

m nd = ' Im ^ ' 
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is the contribution from the M nearest neighbor ions. The total magnetic 
hyperfine field at the impurity site is 

B hf = A(Z imp )m , (14) 

where A(Z imp ) is the Fermi-Segre contact coupling parameter. Similarly we 
define B$ d and 



2. Results 

Our calculations are restricted to compounds where the rare earth ions 
are trivalent. This precludes CeM2, Ybil^and EuM2 compounds: Ce is in 
a mixed valence state and Yb and Eu are divalent. In order to calculate 
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Figure 2: Calculated total magnetic moment mo at Sn impurity diluted in (a) RFe2 and 
(b) RC02 (solid lines) . The dotted lines correspond to the contribution from the rare earth 
and the dashed lines correspond to the contribution of the M ions. 

the local moments and the magnetic hyperfine fields at a Sn or Cd impurity 
diluted in RM2 we have to fix some model parameters. Here, we adopt a 
standard paramagnetic s-p density of states extracted from first-principles 
calculations. The exchange splitting in the s-p energy bands induced by the 
local moments of the rare-earth ions, was properly chosen to yield the s-p 
magnetic moment at the R sites of the host, which is assumed to be of the 
order of 0.1 of the d magnetization at the R sites. The parameter a which 
renormalizes the hopping energy, was chosen a ~ 1 giving the ratio between 
the extension of the host and impurity s-p wave functions. For the whole rare 
earth series we adopted J sd = 0.2 x 10~ 3 in units of the s-p bandwidth in the 
case of -RFe2, while for RC02, we adopted J sd = 0.4 x 10~ 3 for heavy rare- 
earth ions and J sd = 0.5 x 10 -3 for light rare-earth ions both also in units of 
the s-p bandwidth. These choices are enough to reproduce the systematics 
of the respective series. The magnetic moments at a Fe site in _RFe2 ((5' Fe )) 
and at a Co site in RC02 ((S* 00 )) were obtained from Ref. jij. In the present 
case, since we are discussing Laves phase compounds, Z nn = 12. Keeping 
fixed these parameters, we self-consistently determined the local magnetic 
moment and the corresponding magnetic hyperfine field at the Cd or Sn 
impurity. 

Fig. [T] exhibits a typical local density of states (a) at the Cd impurity, 
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and (b) at the Sn impurity diluted at i?Co 2 , namely for R = Gd (Cd is 
not soluble [li] in _RFe 2 ). The local charge potential Vq and the exchange 
interaction with the M ions are such as to produce a slight deformation and a 
shift of the up and down sub-bands generating a total local magnetic moment 
in down direction for the Cd impurity and up direction for Sn impurity (for 
the latter see also Fig. [5] (b)). In fact the Cd impurity local DOS- here 
AZ = 1.2- is close to the rare earth s-p DOS, which has also a negative s-p 
magnetization. On the other hand for the Sn impurity the local potential- 
now AZ = 3.2 -, piles up the down and the up spin states in such way that 
the contribution of the up states to the DOS overcomes that of the down, 
thus generating a change of sign of the magnetic moment with respect to the 
one from the s-p rare earth. 

In table [1] and El the calculated contribution to the magnetic hyperfme 
fields as well the experimental measurements collected from Ref. Q for a Sn 
impurity diluted at R site of the i?Fe 2 and _RCo 2 respectively are displayed. 
The results of a Cd impurity at RCo 2 are shown in the table [31 In both 
tables the agreement with the experimental data are fairly good. 



Table 1: Calculated contributions (from R, induced and total) to the magnetic hypcrfinc 

fields, in Tesla, at Sn impurity in i?Fe2. Experimental data were collected from Ref. Q. 
The mean-square errors is about 0.6 T. 

R Gel Tb Dy Ho Er 

Bf f = 4~5 -3.4 -2.6 -2.3 -2.1 

Bff 56.9 53.7 52.5 49.6 49.3 

B hf 52.4 50.3 49.9 47.3 47.2 

Exp. B hf 52.6 51.2 48.8 47.5 48.0 



Table 2: Calculated contributions (from R, induced and total) to the magnetic hyperfme 
fields, in Tesla, at Sn impurity in i?Co2. Experimental data were collected from Ref. 0]. 
The mean-square errors is about 0.6 T. 

R Gd Tb Dy Ho Er 

Eg -5.3 -2.1 -0.9 -0.4 -0.4 

Bffi 57.2 50.9 45.7 43.9 42.7 

Bhf 51.9 48.8 44.8 43.5 42.3 

Exp. B hf 52.1 47.9 45.0 43.4 41.3 
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Table 3: Calculated contributions (from R, induced and total) to the magnetic hyperfinc 
fields, in Tesla, at Cd impurity in i?Co2. Experimental data were collected from Ref. jfj]. 
The mean-square errors is about 0.09 T. 
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Figure 3: Calculated ratio B^f/n^ at Sn impurity in J?Co2 (full line) and i?Fe2 (dashed 
line). The open circles and full circles represent their respective experimental data dis- 
played in tables [T] and [2] The ratio is almost constant in i?Fe2 but varies in RC02 along 
the rare earth series (see text). 

Figure [2] shows the self-consistently calculated magnetic moments at a Sn 
impurity site in (a) -RFe2 and (b) -RC02. From this figure one can see that the 
contribution to the local magnetic moment arising from the 3d neighboring 
ions, fh l Q d , is larger than the one originated from the rare-earth ions m^, in 
modulus. 

Figure [3] exhibits the _B hf //i 3d ratio, in T//i B , where fj, 3 d is the magnetic 
moment of the 3c? (Fe or Co) atom. We can observe the different behavior 
between Sn:i?Co2 and Sn:i?Fe2 systematics. While in the _RFe2 series the 
Bhi/nzd ratio remains almost constant, in RC02 the ratio decreases from R 
= Gd to R = Er. In both compounds the dominant contribution is from 
the 3d ions (see tables [1] and |2]), the contribution from the rare-earth ions 
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Figure 4: Calculated contributions to the total magnetic hypcrfinc fields Bm at Cd im- 
purity in _RCo2 (solid line). The dotted line correspond to the contribution from the rare 
earth and the dashed from the Co ions, both with the same sign. The squares represent 
experimental data, based on table [3] 
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Figure 5: Calculated magnetic hypcrfinc fields at Cd impurity diluted in i?Ni2. Squares 
represent experimental data. Now the hypcrfinc field varies along the rare earth series. 
See text. 
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being very small. However, the amplitude of the total hyperfine field is about 
twice larger in RCo 2 than in _RFe2 as one goes from Gd to Er. In addition, 
Co has /isd ~ along the series while Fe has about twice larger // 3( j. 
The combination of these aspects, and not an absent contribution from the 
rare-earth sublattice jij, reduce the ratio amplitude in _RFe2 by a factor « 4. 

Fig. H] exhibits the contributions to the total magnetic hyperfine fields at 
the Cd impurity diluted in RCo2- In this case, B^ and B 1 ^ have the same 
sign in agreement with the previous discussion about figure [TJ 

Now, a few remarks concerning the magnetic hyperfine field at Cd impu- 
rity in -RM2, where some experimental data are available Q. In this case, 
the d-band associated to the Ni sublattice is completely filled so no induced 
magnetic moment produced by Ni ions at the Cd impurity is present. Then, 
as illustrated in Figure [5j only contribution due to the rare earth sublattice 
is present. The calculated B^ is thus smaller, in absolute value, than those 
of RCo 2 . 
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